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Abstract. Using new data from the UVES spectrograph on the ESO Very Large Telescope and archive images 
from the Hubble Space Telescope, we have measured projected velocity dispersions and structural parameters for 
two bright young star clusters in the nearby spiral galaxy NGC 5236. One cluster is located near the nuclear 
starburst of NGC 5236, at a projected distance of 440 pc from the centre, while the other is located in the disk 
of the galaxy at a projected galactocentric distance of 2.3 kpc. We estimate virial masses for the two clusters 
of (4.2 ± 0.7) X 10^ Mq and (5.2 ± 0.8) x W'Mq and ages (from broad-band photometry) of lo'^ '^±o.2 ^^^^^ 
"l^gS.oio.i ygg^j-g^ respectively. Comparing the observed mass-to-light (M/L) ratios with simple stellar population 
models, we find that the data for both clusters are consistent with a Kroupa-type stellar mass function (MF). In 
particular, we rule out any MF with a significantly lower M/L ratio than the Kroupa MF, such as a Salpeter-like 
MF truncated at a mass of IMq or higher. These clusters provide a good illustration of the fact that massive, 
globular cluster-like objects ("super star clusters") can form at the present epoch even in the disks of seemingly 
normal, undisturbed spiral galaxies. 
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1. Introduction 

Luminous young stellar clusters have been identified in 
a wide variety of external galaxies, including starburst 
galaxies (van den Bergh 119711 O'Connell et al. 1995J 
Meurer et al. ll995;) . dwarf irregulars (Arp & Sandage .19831 
Melnick et al. 119851 Billett et al. l2UU2fl . merger galaxies 
(Schweizer 120021 Whitmore '20031 , nuclear and circumnu- 
clear starbursts (Harris et al. .20011 Maoz et al. .2001) 
and even in the disks of some normal spirals (Larsen & 
Richtler 11999(1 . Many of these clusters are several magni- 
tudes brighter than any young open cluster known in the 
Milky Way, and the inferred masses can be quite simi- 
lar to those typical of the old globular clusters observed 
around virtually all large galaxies (10"* — lO^Af©). This has 
prompted wide-spread anticipation that observations of 
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such clusters may provide clues to the mechanisms that led 
to the formation of globular clusters in the early Universe. 

One question which has prompted some debate is 
whether or not all of these young luminous clusters will 
actually be able to survive for time spans comparable to a 
Hubble time. While there is ample evidence that at least 
some of them do survive for up to several Gyrs, as il- 
lustrated e.g. by the wide range of cluster ages observed 
in the Large Magellanic Cloud, there have been sugges- 
tions that some young clusters might have a top-heavy 
stellar mass function (MF). Note that we use the term 
'MF' rather than 'IMF' to emphasize that we are observ- 
ing a present-day mass function which may be different 
from the initial one. If a significant fraction of the mass 
is contained in relatively massive stars, the cluster might 
dissolve once these stars reach the endpoint of stellar evo- 
lution (Goodwin ll997|l . However, the most massive young 
star clusters are quite rare, and therefore tend to be lo- 
cated at distances where direct observations of individ- 
ual low-mass stars are impossible with current techniques. 
Thus, the MF can only be constrained by indirect means. 

One way to constrain the MF of stellar clusters is to 
compare a dynamical estimate of the cluster mass with 
theoretical predictions for various MFs, based on simple 
stellar population (SSP) models. Dynamical mass esti- 
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Fig. 1. A VLT/F0RS2 V^-band image of NGC 5236, show- 
ing the locations of the two clusters, and inserts showing 
HST close-ups. Each insert measures 6" x 6". 



mates are usually obtained by measuring the projected ve- 
locity dispersion from an integrated high-dispersion spec- 
trum of the cluster and combining this with information 
about the physical cluster size (typically derived from an 
HST image), assuming virial equilibrium. This method 
was first applied to clusters in NGC 1569 and NGC 1705 
(Ho & FilippenkoCnnSillHniB, and has since been used 
by several other authors in attempts to constrain the MF 
in young clusters. The results so far, however, are far 
from clear. Smith & Gallagher H2001|l found evidence for 
a top-heavy MF in the cluster M82-F, while Maraston 
et al. ((2004) reported "excellent agreement" between the 
dynamical mass and SSP model estimates assuming a 
Salpeter MF extending down to O.IM© for the exceedingly 
massive cluster W3 in the merger remnant NGC 7252. For 
a luminous young cluster in the disk of the nearby spiral 
NGC 6946, Larsen et al. p¥)T)l found a M/L ratio sim- 
ilar to or even greater than expected for a Salpeter-type 
MF extending down to 0.1 M©. Gilbert & Graham 12(^ 
found M/L ratios consistent with "normal" MFs for three 
clusters in NGC 1569. Other authors have found a mixture 
of mass-to- light ratios implying normal MFs in some cases 
and top-heavy MFs in others fSternberg 119981 based on 
a re-analysis of the data by Ho & Filippenko; Mengel et 
al. 12003 McCrady et al. 12003). Evidence for a top-heavy 
MF has also been claimed based on Balmer line strengths 
measured on medium-resolution spectra for clusters in the 
peculiar galaxy NGC 1275 (Brodie et al. lTM^ . 

In this paper we analyse two massive young clusters in 
the nearby spiral galaxy NGC 5236 (M83). In our ground- 
based survey of young massive clusters in nearby spirals 
(Larsen & Richtler HOOl 120001) , NGC 5236 was noted as 
having a very rich population of such clusters, most likely 
because of its high star formation rate. For consistency 
with our previous work, we refer to the two clusters as 
N5236-502 and N5236-805. One of the clusters, N5236- 



Table 1. Template stars. The h values are the peaks of 
the cross-correlation function for the template star spec- 
tra vs. each of the two clusters, using the range 7310- 
7570A. Absolute My magnitudes were calculated using 
apparent V magnitudes from the Bright Star Catalogue 
and Hipparcos parallaxes. 



Star 


Sp type 


h (502) 


h (805) 


Mv 


HR 4226 


Ml II 


0.47 


0.42 


-3.4 ± 1.1 


HR 4279 


Kl II 


0.43 


0.36 


-3.2 ±0.8 


HR 5293 


K4 II 


0.45 


0.40 




HR 5547 


G8 II 


0.43 


0.38 


-1.2 ± 0.3 


HR 5645 


K4 lb 


0.47 


0.42 




HR 5738 


G2 II 


0.41 


0.37 


-2.2 ± 0.7 


HR 6120 


G8 lb 


0.46 


0.41 


-2.3 ±0.8 


HR 6693 


Ml lb 


0.43 


0.42 




HR 7139 


M4 II 


0.42 


0.43 


-2.9 ±0.3 



805, was already noted in our ground-based data, but the 
other (N5236-502) was too red to pass the colour cuts and 
was noted during later inspection of archive images taken 
with the Hubble Space Telescope (HST). The two clusters 
both appear well resolved on HST images, so that their 
half-light radii can be reliably measured without being too 
dominated by resolution effects. We adopt a distance mod- 
ulus for M83 of (m - M)o = 28.25 ± 0.15 (4.5 ± 0.3 Mpc) 
from Thim et al. H2003fl . corresponding to a pixel scale 
on the WF chips of the Wide Field Planetary Camera 2 
on board HST of 2.2 pc pixel^^. Thus, corrections for the 
WFPC2 point spread function (PSF) are still necessary, 
but should not constitute a major uncertainty in the size 
determinations. The HST images show the clusters to be 
reasonably isolated, so we can be confident that the ve- 
locity dispersions measured from the ground are indeed 
those of the target clusters themselves. For measurements 
of the velocity dispersions we use new data from the UVES 
echelle spectrograph on the ESO Very Large Telescope 
(Sec. 121. To determine the cluster ages we use broad-band 
photometry from the Danish 1.54 m telescope at ESO, 
La Silla, supplemented with archive imaging data from 
the FORSl and F0RS2 instruments on the VLT (Sec.|i|). 
The locations of the two clusters within NGC 5236 are 
indicated in Fig. ^ 

2. Velocity dispersions 

2.1. Data 

For the velocity dispersion measurements we obtained 
spectroscopic observations in service mode with the UVES 
echelle spectrograph on the VLT / UT2 (Kueyen) at ESO, 
Paranal, in March 2003. We used a sht width of 0'.'6, 
providing a spectral resolution A/AA ^ 60,000, and a 
dichroic beam splitter to obtain simultaneous exposures 
with the blue side (3730A- 4990A) and red side (66OOA- 
IO6OOA) of UVES. The UVES blue side is covered by a 
single CCD (UVES-B), while the red side uses a mosaic 
of two CCDs (UVES-RL and UVES-RU), resulting in a 



S. S. Larsen and T. Richtler: Dynamical Mass Estimates for Stellar Clusters in M83 3 

Table 2. Velocity dispersions (km/s). Mean values: (1) Weighted by (1/rms), (2) weighted by (h). S/N is per resolution 
element (3 pixels). Wavelength ranges: 4520-4850A (UVES-B), 6770-6850A and 6940-7110A (UVES-RL(l)), 7310- 
7570A (UVES-RL(2)), 7710-8070A (UVES-RL(3)), and 868O-888OA (UVES-RU). 



N5236-502 N5236-805 





A/Tpan 


IVlcLlicLii 








A/Tpan 
iVicdli 


A Aori 1 cm 
IVicU-lcLli 






S /N 




km/s 


km/s 


km/s 






km/s 


km/s 


km/s 






UVES-B 


5.72 


5.77 


0.86 


0.17 


15 


8.46 


8.46 


0.47 


0.12 


27 


UVES-RL(l) 


5.23 


5.15 


0.54 


0.28 


35 


9.57 


9.52 


0.87 


0.17 


36 


UVES-RL(2) 


5.51 


5.42 


0.40 


0.43 


32 


7.82 


7.81 


0.34 


0.39 


34 


UVES-RL(1,2) 


5.36 


5.34 


0.39 


0.35 




8.22 


8.26 


0.31 


0.28 




UVES-RL(3) 


4.90 


5.04 


0.61 


0.28 


27 


7.67 


7.69 


0.59 


0.20 


29 


UVES-RU 


6.33 


6.28 


0.81 


0.34 


24 


7.90 


7.98 


0.82 


0.42 


28 


Mean^ 


5.42 ±0.15 


5.39 ±0.13 








8.13 ±0.18 


8.14 ±0.18 








Mean^ 


5.53 ±0.20 


5.50 ±0.18 








8.02 ±0.19 


8.05 ±0.18 








Adopted 


5.5 ±0.2 










8.1 ± 0.2 












-0.2 r , I , , . I . . . I . . . I . . . I 

-40 -20 20 40 

Pixel shift 



Fig. 2. Illustration of the cross-correlation technique. The 
solid curve shows the cross-correlation function (CCF) for 
cluster N5236-502 vs. the template star HR 4226. The 
three dashed curves show the CCFs for template star HR 
5645 broadened by km/s, 4 km/s and 8 km/s, vs. HR 
4226. For clarity, all CCFs have been shifted to a mean of 
and normalised to a peak value of 1. 

small gap in the wavelength coverage. In addition to the 
two clusters in NGC 5236, one cluster in NGC 2997 was 
observed, but this spectrum turned out to have too low 
S/N for our analysis. We also observed a number of late- 
type bright giants and supergiant template stars, listed 
in Tabled Estimates of the absolute My magnitudes are 
listed for each star, based on the apparent V magnitudes 
in the Bright Star Catalogue (Hoffleit & Jaschek J,99_l ) 
and Hipparcos parallaxes (Ferryman et al. I1997|l . assum- 
ing zero reddening towards these stars. While a parallax is 
available for only one luminosity class lb star (HR 6120), 
it is curious that this star appears to be no more luminous 
than the luminosity class II stars in the table. Part of the 
reason may be that the relative uncertainties on the paral- 
laxes are substantial, as reflected in the large uncertainties 
on the Mv values. Another factor may be our neglect of 
reddening. The Bright Star Catalogue lists a. B — V colour 



of B — V= 1.48 for HR 6120, while G8 supergiants have an 
intrinsic {B-V)o « 1.2 (e.g. Schmidt-Kaler HHHSl- This 
suggests a foreground reddening of E{B — V) ~ 0.3 or 
Ay « 1 mag. The Galactic extinction maps of Schlegel 
et al. H1998|l indicate an extinction of E{B - V) = 0.452 
mag in the direction of HR 6120, which may be considered 
an upper limit to the extinction between us and the star 
itself but is again consistent with Ay « 1 mag towards 
HR 6120. 

The total exposure time for N5236-502 was 11680 s 
(195 min), divided into 4 separate integrations. For N5236- 
805, the exposure time was 7500 s (125 min), divided into 
3 integrations. The initial reductions were done with the 
standard UVES pipeline which bias-subtracts and flat- 
fields the data, extracts the spectra, performs wavelength- 
and flux calibration and merges the echelle orders to a 
single spectrum for each CCD. The final reduced cluster 
spectra were co-added, using a sigma-clipping algorithm 
to remove any remaining cosmic ray hits. Typical S/N val- 
ues for the co-added spectra were ~ 25-35 per resolution 
element (Table EJ. 

2.2. Analysis 

To measure the velocity dispersions of the clusters we em- 
ployed the cross-correlation technique described by Tonry 
& Davis H1979|l . This method has previously been used 
in several other studies of young clusters, including Ho 
& Filippenko ( 1996a! I1996b|l and Larsen et al. plHT]! . 
For the cross-correlation we used the FXCOR task in the 
RV package in IRAF^. Although the UVES pipeline pro- 
duces flux-calibrated spectra, the signature of the indi- 
vidual echelle orders was still clearly visible. To prevent 
these variations from producing spurious features in the 
cross-correlation function (CCF), we fitted a high-order 
spline to the continuum level and filtered out low spatial 

^ IRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc. under contract 
with the National Science Foundation 
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frequencies with a ramp filter before the cross-correlation 
analysis. 

Figure |21 illustrates the basic idea behind our imple- 
mentation of the cross-correlation method: first, the clus- 
ter spectrum is cross-correlated with the spectrum of a 
suitable template star. The fuU-width-at-half-maximum 
(FWHM) of the peak of the resulting CCF is related to the 
broadening of the spectral lines, here assumed to be due to 
the motions of the stars within the cluster. The relation 
between the FWHM of the CCF peak and the velocity 
dispersion is established by broadening another template 
star spectrum with Gaussians corresponding to different 
velocity dispersions, and cross-correlating the broadened 
spectra with the first template. This method assumes that 
the projected velocity distributions are indeed Gaussian, 
which is expected to be a good approximation for re- 
laxed clusters with isotropic velocity distributions (King 
I1965|l . In Fig. 121 the solid curve shows the CCF for clus- 
ter N5236-502 vs. the template star HR 4226. The three 
dashed curves show the CCFs for the spectrum of another 
template star, HR 5645, broadened by km/s (i.e. no 
broadening), 4 km/s and 8 km/s, vs. HR 4226. From this 
example, the cluster appears to have a velocity dispersion 
greater than 4 km/s, but less than 8 km/s. An important 
feature of the cross-correlation technique is that the veloc- 
ity dispersion of the science object is measured relative to 
a template. Therefore, resolution effects and any intrinsic 
broadening of the lines (e.g. due to macroturbulent mo- 
tions in the template star atmosphere) will cancel out, as 
long as these effects are reproducible. Also, because only 
the width of the CCF peak is used in the derivation of 
velocity dispersions, addition of a smooth continuum to 
the spectra (e.g. from early-type stars in the clusters) will 
not affect the measurements. 

The amplitudes of the cross-correlation peaks are 
listed in Table ^ for each template vs. the two clusters, 
cross-correlating over the wavelength range 7310A-7570A 
(UVES-RL(2), see below). The template star spectra were 
broadened with a series of Gaussians with dispersions sep- 
arated by steps of 1 km/s. Except for the three G-type 
supergiants (see below), we used all combinations of the 
template stars in Tabled to produce estimates of the ve- 
locity dispersions. For each combination of template stars, 
the velocity dispersions of the cluster spectra were de- 
termined by interpolation between the dispersions of the 
two best-matching broadened template spectra. Including 
the three G-type stars in the cross-correlation analysis 
significantly increased the scatter, and since the cross- 
correlation signal also tends to be weaker than average 
for these stars we have excluded them from the template 
sample. Presumably the poorer match by the G stars in- 
dicates that the cluster spectra are indeed dominated by 
later-type stars. 

The mean and median values of the velocity disper- 
sions derived for all combinations of template stars are 
listed in Tabled along with the dispersions of the mea- 
surements, the mean peak values of the CCFs (h) and the 
S/N of the cluster spectra. The analysis was carried out for 



several different wavelength ranges: 4520-4850 A (UVES- 
B), 6770-6850 A and 6940-7110 A (UVES-RL(l)), 7310- 
7570 A (UVES-RL(2)), 7710-8070 A (UVES-RL(3)) and 
868O-888OA (UVES-RU). In the UVES-RL(1,2) entry, the 
regions covered by UVES-RL(l) and UVES-RL(2) were 
both included at the same time in the cross-correlation 
analysis. The choice of these wavelength ranges was dic- 
tated partly by the need for adequate S/N (on the blue 
side) and partly by the requirement to avoid spectral re- 
gions that are strongly affected by night sky lines and 
atmospheric absorption bands (on the red side). While 
the UVES blue spectra extend well shortwards of 4500A, 
the contribution of red supergiants to the integrated light 
drops off rapidly at shorter wavelengths. We made some 
crude estimates of the relative contributions of cool and 
hot stars at various wavelengths using Girardi et al. H2()0()|l 
isochrones. The exact ratio is very sensitive to the age 
(and also quite model-dependent), especially for N5236- 
805, but if we define icooi/^tot as the relative fraction of 
the total luminosity contributed by stars with T^ff < 5000 
K, then we find that this ratio decreases from > 50% in 
the /-band to - 25% in V, 5%-10% in the B-band, and 
< 1% in [/ at an age of about 10^ years. Note, however, 
that the cross-correlation analysis does not rely on any 
individual strong lines, but instead utilizes the numerous 
weaker lines present in the spectra of red supergiants. 

2.3. Uncertainties 

It is gratifying to see that the cross-correlation technique 
produces fairly consistent results in the different wave- 
length ranges, with a scatter of less than 10% between the 
different measurements. Even the UVES-B setting pro- 
duces consistent results, although the cross-correlation sig- 
nal here is weak. Furthermore, the scatter of the individual 
measurements within one wavelength range is always less 
than 1 km/s, suggesting that the choice of template star is 
not critical. This is an important point, because the tur- 
bulent velocities in the atmospheres of bright red giants 
and supergiants are comparable (« 10 km/s) to the ve- 
locity dispersions we are trying to measure here (Smith 
& Dominy 1979, Gray & Toner [TWji . The small scatter 
in the velocity dispersion measurements suggests that the 
intrinsic scatter in the macroturbulent velocities among 
the template stars must be relatively small as well (< 1 
km/s). Still, the macroturbulence is an important source 
of uncertainty to keep in mind, and it cannot be excluded 
that the red supergiants in the clusters have systemat- 
ically different properties in this respect than the local 
Galactic stars used as templates. If this is the case, our ve- 
locity dispersion measurements might be biased in either 
direction. For N5236-502, this concern is compounded by 
the fact that the velocity dispersion for this cluster is un- 
comfortably close to the intrinsic broadening of the lines 
in the supergiant stars. For N5236-805 there is another 
problem, namely, that the low age makes it hard to find 
suitable template stars. Stellar isochrones from Girardi et 
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al. H2Ul)U|l predict that the red supergiants have My ~ —6 
at log(age)f»7.1, much brighter than any of the templates 
in Tabled To test whether any strong dependence on the 
template star luminosity class exists, we repeated the anal- 
ysis for the luminosity class lb and II templates separately 
and found that the results agreed within 0.1-0.2 km/s, 
well within the random uncertainty of the measurements. 
We note, however, that a similar analysis of Keck/HIRES 
spectroscopy of four massive clusters studied by Larsen, 
Brodie & Hunter l|20()4(l yielded a somewhat greater dif- 
ference of ^ 0.6 km/s between velocity dispersions de- 
rived from luminosity class la/Ib and II template stars, 
with la/Ib templates yielding lower velocity dispersions. 
That paper also includes a more detailed analysis of the 
dependence of the derived velocity dispersion on the tem- 
plate star properties, including spectral type and luminos- 
ity class. If there is a general trend for the macroturbulent 
velocities to increase with luminosity then we could have 
overestimated the velocity dispersion (and, consequently, 
the mass-to-light ratio) for N5236-805. A more detailed 
investigation of this issue would clearly be of great value. 

As our final estimates of the line-of-sight velocity dis- 
persions we adopted a weighted average of the mean val- 
ues in each wavelength range. We calculated two weighted 
averages for each cluster, one using the r.m.s. dispersion 
within each wavelength range for the weights, and another 
one using the mean amplitude of the cross-correlation 
peak. The result did not depend strongly on the choice 
of weighting scheme and we simply adopted the mean of 
the two weighted averages. For the error estimates, we 
adopted the larger of the estimated standard errors on 
the weighted averages. The final adopted line-of-sight ve- 
locity dispersions are 5.5 ± 0.2 km/s for N5236-502 and 
8.1 ± 0.2 km/s for N5236-805. 

3. Structural parameters 

Structural parameters for the clusters were measured on 
HST images obtained under programmes 5971 (for N5236- 
502) and 8234 (for N5236-814). The data were downloaded 
from the archive at STScI and processed on-the-fly by 
the standard pipeline. The individual exposures were then 
combined with the CRREJ task in the STSDAS package 
in IRAF. For programme 5971, four exposures in F814W 
(2 X 1000s + 2 X 1300s) and two exposures in F606W 
(1100 s + 1200 s) were available, but the central pixels 
of N5236-502 were saturated on all but the short (1000 s) 
F814W exposures, so these were the only exposures used 
for the analysis of structural parameters for this cluster. 
The dataset for programme 8234 consisted of 3 exposures 
in F547M (180 + 350 + 400 s) and 3 exposures in F814W 
(160-1-200-1-350 s), allowing us to carry out separate analy- 
ses in the two bands. Both clusters are imaged on the WF 
chips of their respective HST datasets. 

Interestingly, the HST and VLT images reveal a 
smaller companion cluster, separated by 1'.'2 from N5236- 
502, which is only hinted at in the images from the Danish 
1.54 m telescope as a slight elongation of the cluster im- 
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Table 3. Structural parameters from ISHAPE fits to 
HST/WFPC2 images. Numbers in this table are for a fit- 
ting radius of 3" (30 pixels). For N5236-502 and N5236-805 
the 'V entries refer to F606W and F547M data, while T 
refers to F814W for both clusters. 





N523fi-,'in2 


N.^23fi-Snfi 


-7—. 

Iving models [v) 






-T VV XTiVi 




U. \jO'-±. 






7c _|_ on 








ivlliti iliULieio I J 1 






-T VV nivi 


fY'lzL 

U. i-^ 


U. 




lUU ZL ZU 


Q1 _i_ on 

zn zu 






U. -LOO ZIZ U. UZr 


EFF models (V) 






FWHM 




0'.'084 


Slope parameter 




1.32 


J'off(j"max= OO) 




0'.'140 


reff(rmax= 3") 




0'.'124 


EFF models (7) 






FWHM 




0'.'090 


Slope parameter 


1.14 


1.28 


reff(rmax= Oo) 


I'.'lOl 


0'.'176 


roff(»"max= 3") 


c.'ssg 


0'.'144 


Adopted 






fcff 


0'.'35 ± C.'OS 


0'.'13±0'.'02 



age. The companion cluster is clearly resolved on the HST 
images. The projected separation between the two clus- 
ters corresponds to only 26 pc, so they might constitute 
a real physical pair. From photometry in an r = 5 pix- 
els aperture on the HST images we find a magnitude dif- 
ference between the two clusters of Apgogw = 1-34 mag 
and Af814w = 1-62 mag. Thus, the smaller cluster ap- 
pears somewhat bluer, although the magnitude difference 
in F606W may be slightly underestimated due to the sat- 
uration of the brighter cluster. If the clusters are a physi- 
cal pair, the colour difference might be due to differences 
in the reddenings towards the clusters, which is not too 
unlikely considering that their position coincides with a 
conspicuous dust lane (Fig. Below we estimate an age 
for N5236-502 of about 100 Myr (Section |3|), so the clus- 
ters probably did not form out of material associated with 
this dust lane. 

At the adopted distance of M83, the image scale of the 
WF chips corresponds to 2.17 pc pixel" This is compa- 
rable to the typical half-light radii of stellar clusters, so 
a correction for the WFPC2 point spread function (PSF) 
IS necessary. We used the ISHAPE code (Larsen f TM??)l to 
fit analytic King (jil962(l models and "EFF" (Elson et al. 
I1987f) profiles of the form 

P(r)oc(l + (r/r,)2)-" (i) 

to the images. ISHAPE determines the best fitting shape 
parameters by iteratively convolving the analytic model 
with the PSF, adjusting the FWHM and envelope slope 
parameter a (or concentration parameter c = rt /r,. in the 
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case of King models) until the best fit to the observed 
cluster images are obtained. The PSF was modelled using 
version 6.1 of the TinyTim code (Krist & Hook 1997) and 
the fits were carried out using a fitting radius of 3" (30 
pixels). Integration of Eq. ^ yields the following expres- 
sion for the half-light or effective radius, r^s, of the EFF 
models: 



r,s = rcyil/2)— -I (2) 
for a profile extending to infinity and with a > 1, and 

1/2 



1 



1-a 



1 



- 1 



(3) 



for a profile truncated at an outer radius rmax- For a values 
significantly greater than 1, rctf does not depend strongly 
on the adopted rmax, but this convenient behaviour breaks 
down as a approaches or even decreases below unity. 

The results of the profile fits are listed in Table |21 The 
ISHAPE code does not provide formal error estimates on 
the fitted parameters, and many of the uncertainties are 
likely to be systematic rather than random. An estimate 
of the reproducibility of the results can be obtained by re- 
peating the fits several times with different initial guesses 
for the fit parameters and for different fitting radii. Such 
an exercise shows that the FWHM generally reproduces 
within 0.1 pixel (O'.'Ol) and the EFF a parameter is sta- 
ble within a few percent. The King profile concentration 
parameter is uncertain by at least 20%. In the case of 
N5236-502, the King and EFF fits agree reasonably well 
on the FWHM (0'.'14 and O'.'IT), corresponding to 3.0 and 
3.7 pc. The King fits yield an effective radius of roff = 0'.'36 
(7.8 pc) with a formal uncertainty of about 0'.'05 (1 pc), 
mostly due to the error on the concentration parameter. 
For the EFF fits, the effective radius is sensitive to the 
adopted outer radius of the cluster because the slope a 
is only slightly greater than 1. If we adopt the same rmax 
(3") used for the fit itself, the resulting Toff — 0'.'34. For 
J-max = 2" and r^ax = 5", Eq. Q yields rcff = O'.'SO and 
'''cfi = 0'.'39. In the extreme case of a profile extending to 
infinity, Eq. |5J) yields an effective radius of I'.'IO. For a 
smaller fitting radius of 2" the FWHM increases slightly 
(FWHMw O'.'IS), balanced by an increase in a (- 1.21) 
and leading to a slight net decrease (by about 0'.'03 relative 
to the fits within 3") in rcff. For N5236-805, the envelope 
slope is steeper and Teff is therefore less model-dependent. 
The fits all agree on an effective radius close to O'.'IS. Thus, 
we adopt r^s = 0'.'35 (7.6 pc) and 0'.'13 (2.8 pc) for N5236- 
502 and N5236-805 with uncertainties of O'.'OS (1.1 pc) and 
0'.'02 (0.4 pc). Interestingly, the cluster near the centre of 
NGC 5236 is much more compact than the one located 
further out. 

4. Photometry and cluster ages 

Age estimates were obtained from UBVI colours. TableQl 
lists the integrated V magnitudes and U — B, B — V and 



Table 5. Age and reddening estimates for the clusters. 
For N5236-502 we used the VLT/FORS photometry. The 
Ab values listed here are in addition to our adopted fore- 
ground reddening of Ab = 0.284 mag. 





N5236-502 


N5236-805 




Loff faffed 


Ab 




Ab 




7.99 


0.71 


7.42 


0.24 


Z = 0.008 










UBV 


7.90 


1.31 


6.80 


0.20 


UBVI 


8.06 


1.02 


7.16 


1.00 


Z = 0.02 










UBV 


8.01 


1.13 


7.10 


0.56 


UBVI 


8.11 


0.86 


6.82 


1.07 


Z = 0.05 










UBV 


7.86 


1.42 


7.12 


1.41 


UBVI 


8.06 


0.78 


6.66 


1.55 


Adopted 


8.0 ±0.1 


1.0 ±0.2 


7.1 ±0.2 


1.0 ±0.5 



0.5 



0.0 



-0.5 



-1 .5 



8.5 9.0 

♦n5236-502 




5 j;g;n5236-805 



-0.5 



0.0 



0.5 

3-V)o 



Fig. 3. B~V, U ~B two-colour diagram indicating the lo- 
cation of the two clusters and the S'-sequence of Girardi 
et al. I|1995|l . The logarithm of the age is indicated along 
the sequence. The photometry has been corrected for fore- 
ground reddening only. 



V — I colours for each cluster, corrected for a Galactic 
foreground extinction of Ab = 0.284 mag (Schlegel 
et al. T^M', provided by the NASA/IPAC Extragalactic 
Database (NED)). For both clusters, we obtained pho- 
tometry from images taken with the Danish 1.54 meter 
telescope at ESO, La Silla, using a,n r — 4 pixels (1'.'56) 
aperture (see Larsen 1999 for details concerning the pho- 
tometric calibration). In addition to the data from the 
Danish telescope, wc also list photometry for N5236-502 
from VLT F0RS1/F0RS2 commissioning data (Comeron 
12001(1 . These data have the advantage of much better see- 
ing 0'.'6 rather than 1'.'5), but lack an adequate pho- 
tometric calibration. Therefore, the FORS data were tied 
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Table 4. Ground-based and HST photometry for the clusters. A foreground reddening correction oi As = 0.284 mag 
has been applied. 



Cluster R.A. Decl. 

2000.0 


V 


U-B 


B-V 




y(HST) 
r = 2'.'0 


N5236-502 13:37:06.12 -29:53:18.4 

N5236-502 (VLT/FORS) 

N5236-805 13:37:01.82 -29:52:13.1 


17.51 ±0.01 
17.48 ±0.01 
16.60 ± 0.02 


-0.11 ±0.03 
-0.12 ±0.01 
-0.55 ±0.02 


0.50 ± 0.02 
0.39 ±0.01 
0.19 ±0.03 


0.89 ±0.02 
0.73 ±0.01 
0.65 ±0.04 


17.43 ±0.01 
16.83 ±0.01 




6 7 8 9 

Log(age) 



Fig. 4. Comparison of Bruzual & Chariot model colours 
with observed colours for cluster N5236-502. The hori- 
zontal solid and dashed lines indicate the cluster colours 
corrected for foreground reddening only and assuming 
an additional Ab = 1.0 mag, respectively. The hatched 
areas indicate the ages and adopted uncertainties. The 
model colours are shown for Z = 0.008, Z = 0.020 and 
Z = 0.050. 



to the photometric zero-points of the Danish 1.54 ni data 
using a number of isolated stars outside the main body of 
NGC 5236. The better seeing in the FORS images allows 
a more local measurement of the sky background and a 
smaller aperture radius (4 pixels = 0'.'8) for the cluster 
photometry, reducing possible systematic effects on the 
colours from the nearby companion cluster. Unfortunately, 
N5236-805 is saturated in the FORS images, so for this 
cluster we only have the Danish 1.54 m data to rely on. 
Both the Danish 1.54 m dataset and the F0RS1/F0RS2 
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Fig. 5. Same as Fig. 0] except that cluster N5236-805 is 
shown. 

data include Ha imaging, from which it is clear that nei- 
ther cluster has any line emission associated with it. 

Table 01 also includes estimates of the V magnitudes 
from the HST photometry, obtained in an r = 2" aper- 
ture and using the zero-points in the WFPC2 instrument 
handbook (Biretta et al. ll99H|l . While the central pixels are 
saturated for N5236-502, these contain only a small frac- 
tion of the total flux, and the fact that the 1000 s F814W 
exposures do not reach saturation furthermore suggests 
that the F606W exposures (of comparable length) are not 
strongly saturated. Overall, the V magnitudes from dif- 
ferent sources agree within 0.1-0.2 mag, which probably 
provides a better estimate of the true uncertainty due to 
the background determination, aperture size etc, than the 
formal errors due to photon statistics. It should be noted 
that our errors on the HST photometry do not include un- 
certainties on the photometric zero-points, which may be 
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as large as ~ 0.05 mag according to the WFPC2 instru- 
ment handbook. For the remainder of this paper, we adopt 
the HST values for the total V magnitudes of both clus- 
ters, the VLT F0RS1/F0RS2 colours for N5236-502 and 
the colours from the Danish 1.54 m data for N5236-805. 

In Figure |31 we show the photometry for the two clus- 
ters in a, {B — V , U — B) two-colour diagram, together with 
the "S"-sequence defined by Girardi et al. ( 1995j. The S- 
sequence is basically the mean locus of LMC clusters in 
the {B — V, U — B) plane, and defines an age sequence. 
The tick marks along the sequence indicate the optimal 
direction in which to shift clusters onto the sequence, con- 
sidering reddening effects as well as stochastic colour vari- 
ations. While the S-sequence has certain limitations, such 
as being defined only for LMC metallicity and restricted to 
the B — V, U — B two-colour plane, some rough estimates 
of the ages and reddenings of the two clusters can be drawn 
already from a visual inspection of Fig. El The logarithm of 
the ages of N5236-502 and N5236-805 appear to be about 
7.4 (25 Myr) and 8.0 (100 Myr), respectively For N5236- 
805, only a small shift is required to bring the cluster on 
the S-sequence, suggesting a small reddening correction, 
while a larger shift is necessary for N5236-502. The for- 
mal S-sequence age- and reddening estimates for the two 
clusters are listed in Table |S1 

In order to utilize the full information contained in the 
U BVI photometry and assess the possible role of metal- 
licity differences, we must rely on model calculations for 
the evolution of broad-band colours as a function of age. 
In Figs. ^ and [S] we compare Bruzual & Chariot (2000; 
priv. comm.) model calculations for three different metal- 
licities {Z = 0.008, Z = 0.02 and Z = 0.05) with the ob- 
served colours of the two clusters. In each figure, the hor- 
izontal solid lines indicate the observed U — B, B — V and 
V — I colours (corrected for foreground reddening only), 
while the horizontal dashed lines indicate the effect of an 
additional internal reddening in NGC 5236 oi Ab = 1 
mag. For each metallicity. Table El lists the best-fit ages 
and internal As values obtained by minimizing the r.m.s. 
difference between the model and observed colours, us- 
ing both the full set of UBVI, as well as just UBV. For 
N5236-502, the comparison with model colours yields ages 
that are quite similar to those obtained via the S-sequence 
method, albeit with somewhat higher reddening estimates. 
In the case of N5236-805 there are more pronounced dif- 
ferences, and from Fig. |S1 it is clear that the models are 
very metallicity-sensitive in the age range around the best- 
fit value (~ 10'' years) for this cluster. In particular, the 
"red loop" in the broad-band colours around 10'' years 
is smoothed over in the S-sequence. Thus, the age- and 
reddening estimates are more uncertain for this cluster, 
although the absence of line emission associated with the 
cluster argues against an age of much less than ^ 10'' 
years. 

Based on the various age estimates in Table |5| we 
adopt estimates of log(age) = 8.0 ± 0.1 and 7.1 ± 0.2 and 
reddenings of v4s = 1.0±0.2 and As = 1.0±0.5 for N5236- 



502 and N5236-805. These age estimates are indicated by 
the shaded areas in Figs. ^ and [S| 

5. Modelling mass-to-light ratios 

The Bruzual & Chariot models tabulate mass-to-light 
ratios in different passbands for various stellar MFs. 
However, in order to have the freedom to explore the ef- 
fects of varying the MF and compare results for different 
stellar models, we calculated our own M/L ratios using 
isochrones from both the Padua group (Girardi et al. 12000! 
and references therein) and the Geneva group (Lejeune 
& Schaerer I2001|l . In Fig. we compare various calcula- 
tions of the My magnitude per solar mass for these sim- 
ple stellar population models. As a reference, the solid 
line shows the Bruzual & Chariot SSP models (which are 
based on Padua stellar models) for a Salpeter IjlUSSfl MF 
in the mass range O.I-IOOMq and solar metallicity. The 
dotted-dashed line shows our calculations for a Salpeter 
MF with a minimum mass M,nin = O.lAf0, also based 
on Padua isochrones. The overall agreement between our 
calculations and the Bruzual & Chariot models is quite 
satisfactory and certainly adequate for our purpose. The 
minor differences may well be due to different treatment 
of stellar remnants (we do not count any mass above the 
endpoint of stellar evolution at a given age). The short- 
dashed and dotted lines both show results for a Kroupa 
(,2002.) MF with M,ni„ = O.OIMq, but for Padua and 
Geneva isochrones, respectively. The triple-dotted-dashed 
and long-dashed lines represent Salpeter MFs (based on 
Padua models) with lower cut-offs at IMq and O.QIMq. 
For the youngest ages (log(age)<7.2) there are substantial 
differences between the Padua and Geneva models, but 
otherwise the main difference between the various curves 
lies in the MF choice. 

6. Results 

From application of the virial theorem, the velocity dis- 
persion Vm, total mass M and half-mass radius r^ of a 
star cluster are related as 



fSpitzer 119871 p. 11-12). The 3-dimensional half- mass ra- 
dius r/j is related to the observed half-light radius rcs 
roughly as r^ — 1.3rofr, and with the line-of-sight velocity 
dispersion given as = ^v"^ we can express Eq. (0J in 
terms of observable quantities: 

Mvir « 9.75^^. (5) 

The exact value of the constant in front of Eq. Q de- 
pends on the detailed radial profile of the cluster, but the 
value 0.4 should be accurate to ~10% or better for most 
relevant cases fSpitzer ri987|l . We denote the mass derived 
from Eq. ^ M^ir to distinguish it from a "photometric" 
mass estimate (Mphot), based on SSP models. The above 
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Fig. 6. Absolute My magnitudes per solar mass computed for various combinations of stellar mass functions and 
isochrones. The solid line is from SSP models by Bruzual & Chariot (2000), computed for a Salpeter MF from 0.1- 
lOOM©. The long-dashed, dotted-dashed and triple-dotted-dashed lines are for Salpeter MFs with lower mass limits 
of Mmin = O.O1M0,O.1M0 and l.OM©, computed from Padua isochrones, while the short-dashed hue is for a Kroupa 
MF with Mmin = O.O1M0. The dotted line is again for the Kroupa MF, but using isochrones from the Geneva group. 
All models are for solar metallicity. Data for the two clusters are shown with error bars. 

Table 6. Derived physical parameters. Uncertainties on My are mainly due to the uncertainties on Ab- Photometric 
mass estimates are given for a Kroupa MF and for Salpeter MF with lower mass limit of Q.IMq. Central l^-band 
surface brightnesses (/^o.y) m mag arcsec"^ and estimated core densities (po) in Mq pc~'^ are also listed. 





N5236-502 


N5236-805 


Mv 


-11.57±0.15 


-12.17 ± 0.37 


Log (age) 


8.0 ±0.1 


7.1 ±0.2 


Mphot (Kroupa) [Mq] 


(4.49 ± 0.86) X 10" 


(1.93 ± 1.42) X 10^ 


Mphot (Salp) [Mq] 


(6.56 ± 1.26) X 10'^ 


(2.84 ± 2.06) X 10^ 


^•eff [pc] 


7.6 ± 1.1 


2.8 ±0.4 


A*o,v [uiag arcsec"^] 


14.1 ±0.2 


12.0 ±0.5 


po [Mq pc-^] 


(2.8 ± 1.0) X 10^ 


(1.6 ± 1.1) X 10" 


Mvir [Mq] 


(5.15 ±0.83) X 10^ 


(4.16 ±0.67) X 10^ 


Mphot/Mvir (Kroupa) 


0.87 ±0.22 


0.46 ± 0.34 


Mphot/Mvir (Salp) 


1.27 ±0.32 


0.68 ±0.51 



relations are valid for spherically symmetric clusters with 
an isotropic velocity distribution, which arc in virial equi- 
librium. Another implicit assumption is that the cluster 
light is a good tracer of the underlying mass distribution. 

With these above caveats in mind. Table El compares 
the virial mass estimates (from Eq. 1^1 with photometric 



masses for the two clusters. The photometric mass esti- 
mates are based on Padua isochrones and are calculated 
both for Kroupa and Salpeter (Mmin = O.IMq) MFs. The 
uncertainties on the photometric mass estimates are based 
on the estimated uncertainties on the age (affecting the 
M/L ratio) and on the reddening (affecting My)- These 
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uncertainties are treated as uncorrelated although this is, 
strictly speaking, not true. Errors on the virial masses are 
based on the error estimates for the velocity dispersion in 
Table 121 and on the sizes. The uncertainty on the distance 
to NGC 5236 has not been included in the error budget, 
but the comparison of photometric and virial mass esti- 
mates is not very sensitive to the exact distance (M/L 
ratios scale linearly with distance). The estimated ±0.15 
mag error on the distance modulus quoted by Thim et al. 
(|2()()3|l corresponds to an uncertainty of about 7% on the 
ratio of Mphot/-Wvir- For reference, we also list the central 
V-band surface brightness (/io,v) and the estimated core 
density (po) for each cluster, using the analytic expres- 
sions relating these quantities to the mass-to-light ratios 
and sizes, luminosities and slope parameters of the EFF 
profiles given in Larsen, Brodie & Hunter (,2004,1 . 

For N5236-502, the agreement between photometric 
and virial mass estimates is remarkably good. The virial 
mass estimate falls between the Kroupa and Salpeter pho- 
tometric mass estimates, but is compatible with both 
within the la errors. In Fig. El the cluster actually falls 
right on top of the Kroupa curve calculated for Geneva 
isochrones. An MF with a significant deficiency of low- 
mass stars, such as a Salpeter MF truncated at IMq, 
can be ruled out. For N5236-8G5 the difference between 
the Kroupa and virial masses is somewhat larger than for 
N5236-502. The virial mass is, in fact, somewhat greater 
than expected for a Kroupa MF, implying an excess of 
low-mass stars if the numbers are taken at face value. The 
photometric mass estimates are less secure for this clus- 
ter, however, mainly due to the uncertainties in the stellar 
models around this age, and the inherent uncertainties in 
the age- and reddening determination due to the irregular 
behaviour of the colours and the strong metallicity depen- 
dence, even if the models were perfect. Again, a Salpeter 
MF with a cut-off at IMq can be safely ruled out. The as- 
sumption of virial equilibrium may also be less secure for 
this cluster, due to its young age. We may also have un- 
derestimated our errors on the velocity dispersion, cluster 
size, or both. In conclusion, there is probably no significant 
discrepancy between the virial mass estimate for N5236- 
805 and the assumption that the mass distribution of the 
cluster member stars follows a Kroupa MF. 

How reliable are the velocity dispersions derived from 
integrated spectroscopy of the clusters? One concern is 
that the velocity dispersions might be affected by small 
number statistics, i.e. the integrated light might be dom- 
inated by a few luminous supergiants. To address this 
question, we carried out Monte-Carlo simulations of the 
colour-magnitude diagrams of the two clusters, drawing 
stellar masses at random from a Kroupa MF until the to- 
tal masses equaled the virial mass estimates in Table El 
For each "star" in these synthetic clusters, we then looked 
up the corresponding colour and luminosity in the Padua 
isochrones. For N5236-502 and N5236-8G5, we estimate a 
total of 550 and about 100 red supergiants, respectively. 
The smaller number of red supergiants in N5236-805 is 
mostly due to its younger age and correspondingly higher 



main-sequence turn-off mass (about 13 Mq vs. 5 Mq). 
Thus, we do not expect the velocity dispersions or other 
integrated properties of the clusters to be severely affected 
by small number statistics. 

7. Conclusions 

We have measured velocity dispersions, structural param- 
eters and broad-band colours for two stellar clusters in the 
nearby spiral NGC 5236 (M83). Using these data, we have 
compared photometric mass estimates, based on Padua 
stellar isochrones and various assumptions about the stel- 
lar mass function (MF), with virial mass estimates for 
the two clusters. We conclude that the observed mass-to- 
light ratios are consistent with a Kroupa-type MF, and we 
can rule out an MF with a significant deficiency of low- 
mass stars relative to the Kroupa MF. We derive virial 
mass estimates of (5.2 ± 0.8) x IO^Mq (N5236-502) and 
(4.2 ± 0.7) X IO^Mq (N5236-805), comparable to or even 
somewhat higher than the typical masses of old globu- 
lar clusters. The velocity dispersions of these clusters (5.5 
km/s and 8.1 km/s) are probably close to the limit at 
which it becomes impractical to measure them from inte- 
grated light. For even lower velocity dispersions, system- 
atic uncertainties due to variations in the turbulent ve- 
locities of the red supergiant stars, mass segregation, and 
other factors, will pose a formidable challenge. 

While one of the clusters is located near the nucleus of 
NGC 5236 and may disrupt on relatively short timescales, 
the other is located in the disk of the galaxy at a projected 
galactocentric distance of 1'.8, corresponding to 2.3 kpc, 
and has an estimated age of about 100 Myrs. The masses 
of these clusters are comparable to those of the luminous 
young clusters (sometimes referred to as "super star clus- 
ters") encountered in merger galaxies and starbursts, but 
they are much higher than that of any young open cluster 
observed in the Milky Way today. These clusters provide 
a striking example of the fact that such clusters can still 
form in the disks of some apparently normal, undisturbed 
spiral galaxies. 
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